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Much  of  the  work  here  lias  been  computationally  oriented,  focusing  on  a 
variety  of  problems  involving  a  varying  region  whose  boundary  is  to  be  deter¬ 
mined  (i.c  ,  free  boundary  problems  and  shape  optimization).  The  compu¬ 
tational  approaches  have  included  use  of  boundary  tracking  with  remeshing, 
applying  an  auxiliary  function  on  a  fixed  grid,  and  reformulation  as  a  bound¬ 
ary  element  method  (BEM). 

Collaboration  has  continued  with  a  group  at  Universite  de  Paris-Sud  (Or- 
say)  on  theoretical  and  computational  approaches  to  a  problem  of  tracking 
the  evolution  of  a  solid/liquid  interface  at  which  one  has  accretion  or  dissolu¬ 
tion  of  a  substance  dissolved  and  diffusing  (possibly  with  chemical  reactions) 
in  the  liquid  phase.  Both  fixed  and  adaptive  mesh  algorithms  have  been  de¬ 
veloped  and  are  to  be  implemented  for  comparative  study.  Of  interest  are 
the  computational  technique  itself,  the  average  intcrfacial  velocity,  and  the 
stability  of  the  shape. 

A  model  problem  for  structural  ••plimizatiou  is  bring  studied  jointly  with 
C.  Schwab.  After  derivation  of  Du  variation  for  the  relevant  functional  a*  a 
boundary  integral,  this  is  being  i  •  mpulnlionallv  l rented  via  |)EM.  For  the 
simplest  case  (only  an  area  constraint)  the  exact  solution  is  known.  'Hie 
method  has  been  implemented  for  '!ii«  and  rapidly  finds  the  optimal  shape 
for  rather  general  initial  shapes  Two  interesting  difficulties  arise  for  which 
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investigation  continues:  (i)  since  the  computational  boundary  is  inevitably 
polyhedral,  one  must  treat  corner  singularities  and  this  has  proved  surpris¬ 
ingly  tricky  in  practice,  (ii)  maintaining  accuracy  can  be  expected  to  require 
occasional  retneshing  and  so  alteration  of  the  dimensionality  of  the  outer 
optimization  problem.  The  latter  question  also  occurs  in  consideration  of  a 
coordinated  approach  to  refinement  of  discretization  and  optimization,  now 
under  analysis. 

Other  activity  includes  work  on:  (i)  switching  systems,  (ii)  ill -posed  prob¬ 
lems,  (iii)  collaboration  with  S.  Anttnan  on  a  problem  of  nonlinear  viscoHas- 
ticity  (longitudinal  motions  of  a  rod),  and  initiation  of  collaboration  with 
S.  Lcnharl  on  biorcmcdiation.  Activity  in  (i)  includes  a  Ph.l)  degree  for 
W.  Szczcchla,  extension  of  results  by  Filippov,  Utkin,  etc.,  on  sliding  modes 
arising  in  hybrid  (variable  mode)  systems,  and  collaboration  with  M  A  Kras 
nosel’skTi  on  existence  of  periodic  solutions  of  forced  systems.  Work  in  (iii) 
has  now  led  to  a  proof  that  a  certain  inequality  constraint  on  the  constitu¬ 
tive  relation  implies  the  impossibility  of  infinite  compression.  Finally,  work 

in  (iv)  has  resulted  in  an  interesting  formulation  of  a  biorcmcdiation  problem 
involving  ODBs  with  discontinuous  right  hand  side,  as  in  (i);  the  optimal  con¬ 
trol  problem  is  now  under  consideration,  as  is  the  extension  to  a  Distributed 
Parameter  problem  involving  flow  of  ground  water  and  a  free  boundary  prob¬ 
lem. 


Marc  Teboulle 


Within  the  general  framework  of  the  proposal,  work  on  nonlinear  opti¬ 
mization  has  been  primarily  in  the  areas  of:  (i)  dual-based  methods  and 
decomposition,  (ii)  regularization  and  approximation  techniques,  (iii)  non*- 
month  optimization,  and  (iv)  algorilm  development  for  large-scale  optimiza¬ 
tion.  These  methods  arc  applicable  to  o  wide  range  of  important  applications 
including,  e.g.,  optimal  shape  design,  structural  optimization,  and  image  re¬ 
construction.  Since  the  high  dimensionality  of  these  problems  poses  severe 
computational  challenges,  work  lias  focused  on  dual-based  and  augmented 
Lngrangian  methods  taking  advantage  of  available  special  structure  to  permit 
natural  decompositions  for  which  parallel  iterative  methods  can  be  utilized 
elfec  lively. 
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Iii  particular,  this  methodology  has  been  applied  to  problems  of  structural 
optimization  and  image  reconstruction.  Initial  analysis  and  supporting  tech- 
nical  material  are  described  in  a  sequence  of  papers,  partly  in  collaboration 
with  a  doctoral  student  G.  Chen  and  with  A  lusem.  Related  work  on  approx* 
im alien  of  a  class  of  optimization  problems  and  application  of  optimization 
methods  to  problems  of  robust  stabilization  has  been  in  collaboration  with 
.1.  Kogan  and  with  i.  Vajda. 

Papers  by  Tebouile  and  by  Chen  and  Teboulle  introduce  the  basic  theory 
and  convergence  analysis  for  constructing  a  general  new  class  of  augmented 
Lagrangians.  The  approach  is  giving  rise  to  new  multiplier  methods  which 
appear  to  have  many  attractive  characteristic  properties  and  which  arc  use¬ 
ful  for  decomposition;  iti  particular,  these  arc  as  differentiable  as  the  data 
permits,  opening  the  possibility  of  applying  second-order  methods  after  de¬ 
composition  to  smaller  problems. 

A  byproduct  is  a  way  to  generate  highly  parallelizable  methods  for  so¬ 
lution  of  discrete-time  deterministic  control  problems.  Some  resemblance 
to  recent  work  of  Rockcfcllar  suggests  the  possibility  of  extension  to  linear- 
quadratic  problems  arising  in  multistage  stochastic  programs,  etc.  Initial 
results  will  appear  in  a  paper  with  Chen  This  approach  should  be  particu¬ 
larly  useful  for  structural  optimization  problems  which,  as  typically  oerur*, 
have  a  large  but  tractable  inner  problem  and  a  smaller  but  difficult  (mms- 
nioolh)  outer  problem. 

Another  application  is  to  nn  optimization  formulation  of  image  recon¬ 
struction.  In  a  paper  by  Teboulle,  a  convex  regularization  process  for  linear 
inverse  problems  is  treated  via  duality  methods.  An  iterative  method  for  a 
specific  image  reconstruction  from  noisy  projections  wits  studied  with  lusem, 
using  a  regularization  of  the  objective  functional,  enforcing  smoothness  to 
improve  the  quality  of  the  images 


